Controlling the Kondo Effect in CoCu n Clusters Atom by Atom 
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Clusters containing a single magnetic impurity were investigated by scanning tunneling mi- 
croscopy, spectroscopy, and ab initio electronic structure calculations. The Kondo temperature 
of a Co atom embedded in Cu clusters on Cu(lll) exhibits a non-monotonic variation with the 
cluster size. Calculations model the experimental observations and demonstrate the importance of 
the local and anisotropic electronic structure for correlation effects in small clusters. 
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Nanometer scaled electronic devices require the under- 
standing and the control of electron behavior, in partic- 
ular of correlation effects, at the atomic scale. Experi- 
mental techniques such as scanning tunneling microscopy 
(STM) and spectroscopy and angle-resolved photoemis- 
sion demonstrate the relevance of many-body phenomena 
beyond standard band theory [1]. In many cases, espe- 
cially for compounds of d and / elements, strong elec- 
tron correlations should be taken into account even for 
a qualitative description of electron energy spectra and 
physical properties Q. The Kondo effect, which is re- 
lated to the resonant scattering of conduction electrons 
by quantum local centers, is one of the key correlation ef- 
fects in condensed matter physics [3]. It arises from the 
interaction of a single localized magnetic moment with a 
continuum of conduction electrons and results in a sharp 
Abrikosov-Suhl resonance of the local impurity spectral 
function at the Fermi level and below a characteristic 
Kondo temperature, Tk 01 • Originally detected as a re- 
sistance increase below Tk of materials with dilute mag- 
netic impurities, this many-body phenomenon has been 
observed for semiconductor quantum dots 0, H[ and for 
atoms 0, SEilH and molecules at surfaces. Whereas 
the understanding of correlation phenomena of bulk ma- 
terials and single impurities has made significant progress 
much less is known about these effects in structures such 
as clusters where typical dimensions are on an atomic 
scale. STM offers unique opportunities to fabricate such 
structures by manipulation of single atoms and a few at- 
tempts have been made to achieve a degree of control 
over the electronic |U E El and magnetic structure 
of atoms and clusters on surfaces. 

Here we show that the Kondo effect of a magnetic 
impurity does not simply scale with the number of 
nearest neighbors. Rather, the local and anisotropic 
electronic structure at the impurity site is important 
for correlations in atomic scale structures, where "each 
atom counts" [16]. We investigate the evolution of the 
Abrikosov-Suhl resonance induced by a Co atom ad- 
sorbed on Cu(lll) as we change the number of nearest 



neighbor Cu atoms. A non-monotonic behavior of the 
Kondo temperature with the cluster size is at variance 
with bulk models which are typically invoked to interpret 
the Kondo effect of a single adsorbed atom (adatom). 
Scanning tunneling spectroscopy measurements and ab 
initio calculations clearly reveal that the Kondo effect 
in atomic scaled clusters depends crucially on their de- 
tailed geometric structure and thus on the local density 
of conduction electron states at the magnetic site. 

Experiments were performed with a home-made scan- 
ning tunneling microscope operated at 7 K and in ultra- 
high vacuum with a base pressure of 10 -9 Pa. Tung- 
sten tips and Cu(lll) surfaces were prepared by argon 
ion bombardment and annealing. While single Co atoms 
were deposited onto the cold surface using an electron 
beam evaporator and an evaporant of 99.99 % purity, 
single Cu atoms were transferred from the tip as pre- 
viously reported |ljj. The individual atoms were chemi- 
cally identified by the presence (Co) or absence (Cu) of 
the Abrikosov-Suhl resonance. Clusters consisiting of a 
single Co atom and several Cu atoms HI, [Hj were fab- 
ricated by atom manipulation with the microscope tip. 
Spectroscopy was performed by a state-of-the-art lock-in 
technique. 

Figures [TJ a) -(e) show a series of constant-current STM 
images, which illustrates fabrication of a CoCu cluster 
[Figs.QJa), (b)] and dimensions of CoCu n (n = 1, . . . ,4) 
clusters [Figs.QJb)-(e)]. In STM images acquired at a 
sample voltage of 100 mV Co atoms appear higher than 
Cu atoms [Fig.[TJa)], and thus an additional means of 
discriminating the adatom species is provided. The stick- 
and-ball models of the clusters used for the calculations 
are presented in Fig. [Iff). Spectroscopy of the differen- 
tial conductance (dl/dV) performed with the tip posi- 
tioned above the cluster center reveals an unoccupied 
state whose energy decreases with increasing number of 
Cu atoms [squares in Fig.[]Jg)]. Our calculations show 
that this resonance is of p z character and are in agree- 
ment with experimental data [circles in Fig.QJg)]. This 
resonance likewise serves as an indicator of the cluster 
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FIG. 1: (Color online) Scanning tunneling microscopy images 
and resonance energies of CoCu n clusters acquired at 7K. (a) 
STM images of single Co and Cu atoms adsorbed on Cu(lll) 
prior to fabricating CoCu n clusters, (b) - (e) Sequence of 
STM images showing CoCu n (n = 1,...,4) clusters. The 
same lateral and height scale was used for each image. Images 
were acquired at a sample voltage of 100 mV and a tunneling 
current of 0.1 nA. (f) Stick-and-ball models of fully relaxed 
adsorption structures of CoCu n (n = 0,...,4) on Cu(lll). 
(g) Experimental energy of unoccupied cluster state versus 
the cluster size. The size of the symbols (squares) corresponds 
to the uncertainty of the energies. Calculated energies are 
depicted as circles. 



composition and size. 

Figure [2] shows a sequence of dl/dV spectra acquired 
on a single Co atom (lower curve, n = 0) and above 
the center of CoCu n clusters with n ranging between 1 
and 4. Starting from a sharp indentation of the differ- 
ential conductance close to the Fermi level (sample volt- 
age V = OmV) which is the spectroscopic signature of 
the single-Co Kondo effect on Cu(lll), the dl /dV curve 
broadens appreciably upon adding two Cu atoms to the 
Co atom. Surprisingly, upon increasing the number of hy- 
bridizations to three and four Cu atoms, dl/dV spectra 
exhibit a sharpening of the resonance again in contrast 
with the monotonic behaviour of the unoccupied p z -like 
state with the number of copper atoms [Fig.QJg)]. 

To quantify the broadening of the resonances at the 
Fermi level we described experimental data by a Fano 
line shape [Hj]: 



dl (q 



(i) 



with q the asymmetry parameter of the Fano theory and 
e = (eV — ckV^bTk) where — e denotes the electron 
charge, V the sample voltage, ex the resonance energy, 
and ke Boltzmann's constant. The profiles of Fano line 
shapes according to Eq.[T] and fitted to dl/dV data are 
presented as solid lines in Fig.[2j The striking behavior of 
broadening and sharpening of the Abrikosov-Suhl reso- 
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FIG. 2: (Color online) Spectroscopy of the differential con- 
ductance (dl/dV) of CoCun clusters around the Fermi level 
performed at 7K. The spectroscopic feature is the Abrikosov- 
Suhl resonance induced by the Kondo effect. Solid lines depict 
fits of a Fano line shape (Eq.[TJ) to experimental data. Prior 
to spectroscopy the tunneling gap was set at 30 mV and 1 nA. 
Spectra for n — 1,2,3,4 are vertically offset by 5, 10, 20, 
25 nS, respectively. 



TABLE I: Experimental and calculated Kondo temperatures 
(Tk) and asymmetry parameters (q) for clusters CoCu n . Tk 
and q were extracted from fits of the Fano line shape (Eq.[T]) 
to measured Kondo resonances shown in Fig. [2] Experimen- 
tal uncertainty margins correspond to standard deviations ob- 
tained by a statistical analysis of a variety of fits. Calculated 
Tk and q are the arithmetic mean of values obtained by a 
Gaussian broadening of the sp local density by 100 meV and 
50 meV. The uncertainty margins reflect the deviations of the 
arithmetic mean from data obtained by 100 meV and 50 meV 
broadening. 



Experiment 



Calculation 
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T K (K) 


q 


T K (K) 


q 





61 ±4 


0.09 ± 0.03 


52 ±4 


0.12 ±0.01 


1 


142 ± 10 


0.12 ±0.03 


204 ± 12 


0.03 ±0.01 


2 


326 ± 30 


0.07 ±0.03 


288 ±3 


0.01 ±0.01 


3 


200 ± 15 


0.22 ±0.04 


221 ± 18 


0.13 ±0.01 


4 


43 ±6 


0.22 ±0.04 


79 ±3 


0.33 ±0.01 



nance with increasing number of Cu atoms is reflected by 
the Kondo temperature Tk which we plotted in Fig.[3](a). 
While the single Co adatom exhibits Tk = (61 ± 4) K, 
which is in agreement with a previous publication [9[ , for 
C0CU2 we find (326 ± 30) K which then decreases again 
to (43 ± 6)K for CoCu 4 . Figure Mb) shows the evolu- 
tion of the asymmetry parameter with the cluster size. 
For n = 0,1,2 the asymmetry parameter varies weakly 
around 0.1. It increases steeply by ~ 100 % upon adding 
a third Cu atom. Kondo temperatures and asymmetry 
parameters for the various clusters are summarized in 
Table [D 

The peculiar behavior of the Kondo resonance with 
cluster size is at variance with the monotonic dependence 
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FIG. 3: (Color online) Experimental (squares) and calculated 
(circles) Kondo temperatures and asymmetry parameters, (a) 
Kondo temperature (Tk) as a function of cluster size, (b) 
Asymmetry parameter (q) as a function of cluster size. 



of Tk on the average hybridization strength expected 
from bulk models of the Kondo effect. Theoretical ap- 
proaches to the Kondo effect of magnetic impurities on 
surfaces consider the importance of bulk 

[13, and sur- 
face |22|, [23|, [24[ states in the scattering of electron waves 
at the magnetic impurity site. Recent experiments with 
Co adatoms on Cu(100) [9[, which does not host any sur- 
face state close to the Fermi level [25|, or on Cu(lll) 
but in the vicinity of atomic surface steps that affect the 
surface states [26| have indicated the importance of bulk 
rather than surface states for the Kondo effect. In this 
situation, ab initio electronic structure calculations which 
are not restricted by specific model assumptions can be 
very useful. Here, we perform accurate calculations of 
the local electronic structure of clusters on the copper 
surface. 

The Kondo temperature for a single magnetic impurity 
can be estimated as 



T K w W y/\JN(E F )\ exp 



\JN(E ¥ 



(2) 



where W is the conduction (s) electron bandwidth, J the 
s-d exchange integral, and N(E-p) the density of states 
at the Fermi energy Ey [3[. Since the s-d exchange in- 
teraction is local, N(E-p) in Eq.[2]is the local density of 
states of conduction electrons at the magnetic atom. 

The asymmetry parameter, g, is also determined by 
the local electronic structure. In a simple model q may 
be expressed as 
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FIG. 4: (Color online) Calculated local density of states 
(LDOS) at the Co site for CoCu n assemblies (n — 0, . . . , 4). 
(a) sp LDOS plotted for energies between — 1.9eV and 1.9 eV 
with respect to the Fermi energy, (b) Close-up view of (a) 
showing the energy scale of the Kondo resonance, (c) p z 
LDOS with indicated energies of unoccupied cluster reso- 
nances shown in Fig.[TJg). 



where G is the local Green's function of the conduction 
electrons at the impurity site and 7 measures the ratio of 
the coupling of the scanning tunneling microscope tip to 
conduction electron states and to the strongly localized 
Co d states [27j. 

To realistically describe the CoCu n clusters on Cu(lll) 
we calculated their electronic and structural properties 
by means of density functional theory using the gener- 
alized gradient approximation to the exchange correla- 
tion potential [28]. For solving the resulting Kohn-Sham 
equations accurately, the Vienna Ab Initio Simulation 
Package J29[ with the projector augmented wave basis 
sets |30|, ISJj and 350 eV as plane wave cut-off have been 
used. We modeled the CoCu n structures using 4 x 3 su- 
percells of Cu(lll) slabs containing up to 7 Cu layers. 
These structures haven been fully relaxed with the re- 
quirement that all forces are less than 0.2eVnm _1 and 
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were then used to calculate the local density of states. 
The supercell Brillouin zone integrations for obtaining 
the local density of states have been performed using the 
tetrahedron-method on 6 x 6 k meshes with subsequent 
50 - 100 meV Gaussian broadening. Below, we use both, 
50 and 100 meV smearing, to provide an uncertainty mar- 
gin for the derived Kondo temperatures and asymmetry 
parameters. 

The relaxed structures [Fig.QJf)] show that the dis- 
tance of the Co atom to the Cu surface increases mono- 
tonically upon adding Cu atoms to the cluster. However, 
the local electronic structure at the Fermi level and at 
the Co site varies in a non-monotonic way with increas- 
ing number of Cu atoms. The substrate conduction band 
states extend to the Co adatom site where they interact 
with localized Co 3d electrons. This interaction is quanti- 
fied by the Co sp-projected local density of states, N(E), 
which is shown in Figs.[D(a) and[D(b). The dependence 
of G(E-p) on the number of Cu atoms is non-monotonic 
[see the close-up view around in Fig. [4|(b)] and leads, 
according to Eq.0 to the experimentally observed non- 
monotonic behavior of the Kondo temperature (Fig. [3]). 
We used W = 20 eV and J = 1.3 eV as fitting parame- 
ters. Moreover, with 7 = 0.22 eV -1 in Eq.[3]the variation 
of G(E-p) with the number Cu atoms also reproduces the 
non-monotonic trend of q observed in the experiments 
[Fig. [31(b)]. Figure Htc) shows the evolution of unoccu- 
pied resonance energies with the number of Cu atoms. 
The resonances appear as peaks in the calculated p z - 
projected local density of states whose energies show the 
same monotonic behavior as in the experiment [Fig.[]Jg)]. 

The non-monotonic dependence of N(E-p) on the num- 
ber of Cu atoms surrounding the Co atom is related to 
the specific nature of chemical bonds in the system. Since 
localized d orbitals of Co atoms are strongly anisotropic 
and form well-defined directional bonds the degree of 
their hybridization with Cu atoms is determined by the 
whole geometry of the cluster rather than by just the 
coordination number. 

In summary, we have shown that the Kondo effect in 
clusters is but depends crucially on their detailed geo- 
metric structure. Atom by atom manipulation changes 
the local density of conduction electron states at the 
magnetic site and thus varies the Kondo temperature 
in rather broad limits. The results demonstrate that 
electron correlations may be tuned in atomic-scale struc- 
tures. 
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